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Abstract 



We examine the effects of electron-phonon interaction on the dynamics 
of the charge carriers doped in two-dimensional (2D) Heisenberg antiferro- 
magnet. The t-J model Hamiltonian with a Frohlich term which couples the 
holes to a dispersionless (optical) phonon mode is considered for low doping 
concentration. The evolution of the spectral density function, the density of 
states, and the momentum distribution function of the holes with an increase 
of the hole-phonon coupling constant g is studied numerically. As the cou- 
pling to a phonon mode increases the quasiparticle spectral weight decreases 
and a "phonon satellite" feature close to the quasi-particle peak becomes 
more pronounced. Furthermore, strong electron-phonon coupling smears the 
multi-magnon resonances ( "string states" ) in the incoherent part of the spec- 
tral function. The jump in the momentum distribution function at the Fermi 
surface is reduced without changing the hole pocket volume, thereby provid- 
ing a numerical verification of Luttinger theorem for this strongly interacting 
system. The vertex corrections due to electron-phonon interaction are negli- 



gible in spite of the fact that the ratio of the phonon frequency to the effective 
bandwidth is not small. 
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I. INTRODUCTION 



Since the discovery of high-temperature superconductivity in copper-oxide based mate- 
rials, experimental studies have revealed a lot of evidence that interaction between electrons 
and lattice vibrations plays an important role in these compounds.S The changes in position 
and width of the phonon peaks below superconducting transition temperature measured by 
Raman spectroscopy^ and by neutron scattering techniquesi, the existence of the isotope 
effect which varies as a function of both doping and rare earth ion substitutionB, clearly 
demonstrate that coupling between charge carriers and phonon modes in high-T c supercon- 
ductors is not negligible. Nevertheless, theoretical studies of the effects of electron-phonon 
interaction on the dynamics of the charge carriers in electron-correlated systems are far 
from complete. The main problem for proper treatment of electron-phonon interaction in 
high-T c superconductors is that charge carrier motion in an antiferromagnetic background (or 
spin-fluctuating background at higher doping concentration) is strongly affected by electron- 
magnon interactions itself.lH0 Including the electron-phonon interaction deals with quasi- 
particles renormalized by interaction with magnons.0 As a consequence the renormalized 
bandwidth becomes comparable to phonon frequencies and the 'classical' Migdal-Eliashberg 
approach to the electron-phonon problem^ in metallic systems seems beyond the region of 
application. 

In the present paper, we add a Frohlich term to the t — J model Hamiltonian to study 
the quasiparticle properties in the presence of hole-phonon interaction. This is particularly 
interesting, since it is widely believed that the t — J model incorporates essential features 
of electron systems with strong local repulsion between electrons, characteristic of high T c 
copper oxides.0 We examine the properties of this model at low doping concentrations (near 
half-filling) where analytical and numerical results for the t — J model are well-known. The 
optical phonon mode in the copper-oxide plane is assumed to have a substantial influence on 
the dynamics of doped holes in cuprate compounds. In the past, quasiparticle properties in 
the presence of electron-phonon interaction were studied by Engelsberg and SchrieffeJi for 
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a weakly correlated model of conventional metal. The authors examined the spectral density 
function A(k,u) of an electron for phonon spectra of the Einstein and Debye forms. They 
found that spectral density function exhibits several branches of excitations rather than a 
single branch of a dressed electron. Our calculations of the quasiparticle energy dispersion 
show the same for the present model of a strongly correlated system with hole-phonon 
interactional. Recently the mass renormalization due to a coupling to optical phonons was 
studied for strongly correlated electrons by Ramsak et a/.0. The authors found that the 
phonon-induced mass renormalization of a single hole that propagates in the t — J model on 
the scale 2 J is much larger than that in the corresponding uncorrelated model for J <^ t. The 
mass enhancement, increasing with J/t ratio, is due to the slow motion of a spin polaron, 
which makes hole-phonon interaction more effective. The hole excursions from the center 
of the spin polaron are restricted to a few lattice sites and the bandwidth is on scale t 2 /J, 
as long as J/t is large. On the other hand, when J <C t the confining antiferromagnetic 
potential becomes weak and a hole performs large radius incoherent excursions on scale t. 
Therefore, the mass enhancement induced by the hole-phonon interaction becomes weaker 
and the coherent part of the spectral weight tends to zero. 

In our paper, we study various quasiparticle properties as well as the mass renormaliza- 
tion due to hole-phonon interaction in the limit of a finite but small concentration of doped 
holes. The experimentally measured optical phonon frequencies!^ are close to 50 meV and 
typical values of J are of the order of 100 meV . We assume the optical phonon frequency, Q, 
to be Q = 0.5 J. The t — J model Hamiltonian appended with Frohlich term will be denoted 
as the t — J — g model. Our paper is organized as follows. In Sec. [T| the quasiparticle residue 
and the mass renormalization constant are calculated for t — J — g model. The incoherent 
part of the spectrum and Luttinger's theorem are studied in Sec. |TTI|. Vertex corrections 
due to the hole-phonon interaction are analyzed in Sec. [TV| Section [V] presents the optical 
conductivity in the t — J — g model. A summary is given in Sec. VI. 
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II. FORMULATION AND QUASI-PARTICLE SPECTRAL WEIGHT 



To analyze the interaction between optical phonons and charge carriers in the copper 
oxide planes, we consider a two-dimensional t— J model Hamiltonian appended by a Frohlich 
term0 with the dispersionless (optical) phonon mode f2 and the coupling constant g. Using 
linear spin-wave approximation through Holstein-Primakoff transformation^! and spinless- 
fermion/ Schwinger-boson representationi of the electron operators on a site i with spin a 
Cia- cji — > hiS^ and — > h\ the Hamiltonian can be written: 

HtJg = -* E [ h i k J S i + S i h i k t + k A + J E P'A ■ SjPj 

{id) (ij) 

-fi^Pi+aJL KHh + bf) + a E Kh , (i) 

i i i 

where hi is a hole creation operator on site i and Si is the on-site electron spin operator. 
(5j r turns down the spin of the electron on site i.) The summation is over the nearest 
neighbor sites of the square lattice, bf is an optical phonon creation operator at site 
i. p\ = hfhi is the local density operator of the spinless hole, while p\ = 1 — p\ is the 
local density operator of electron under the single occupancy constraint. Further, p is the 
chemical potential related to the concentration of the doped holes: (p%) = x (at half-filling 
x = 0). 

Since t > J would correspond to a situation in cuprates, the usual perturbation ap- 
proximation does not work in the strongly coupled t — J model. However, it has been for 
some time known that the self-consistent noncrossing approximation for the hole self-energy 
gives a fairly accurate result because of the small vertex corrections. Namely, it was ex- 
plained in Refii!, that the leading correction to the hole-spin wave vertex vanishes as a 
consequence of the electron (hole) spin conservation in the t — J model. Mathematically, 
this fact is reflected in the special symmetry of the hole-spin wave interaction vertices M^^ 
and Ngff, which leads to the identical vanishing of non- maximally crossing hole-spin wave 
vertex corrections^. However, it is uncertain whether or not the vertex corrections due to 
the hole-phonon interaction are negligible in a strongly correlated electron system. Thus, 



in the first place we assume that they are so, and later in Sec. [TV] this assumption will 
be supported by calculating the lowest order vertex corrections. To study the influence of 
hole-phonon coupling on the quasiparticle properties we treat the Hamiltonian within the 
noncrossing approximation for spin wave and phonon interactions^. Fig. |1] shows four possi- 
ble noncrossing diagrams contributing to the self-energy: spin wave emitting and absorbing 
diagrams, and phonon emitting and absorbing ones. The thick solid line is the dressed hole 
Green's function, while the unperturbed hole propagator is denoted by the thin solid line. 
The dashed and wavy lines stand for the spin wave and phonon propagators, respectively. 
Here M, N are the hole-spin wave interaction vertices, while g is the hole-phonon vertex. 
After summing over the intermediate frequency, the hole self-energy at zero temperature 
becomes 

A(k-q,y) at2 ro ^ A{k-q,y) 



E R (k,u)=y2{MlJ°°dy A[k q ' V) S +N$J dy 



roc 



+ / / dy * k -™. s + j r dy A{k :f: y) . s } , (2) 



lo uj — y — Q + i5 J-oo uj — y + Q + iS 
where the superscript R indicates retarded functions which are analytic in the upper half- 
plane of u, Eq is the spin wave energy, and A(k,u) is the spectral density function of the 
interacting hole propagator given as 

A(k,u) = --ImG R {k,uj) . (3) 

7T 

The hole Green's function is then obtained from the self-energy 

G R (k,u) = —J- . (4) 

uj — L u (k, uj) + jj + to 

The self-consistent integral equation, Eq. [|, is solved numerically (See Ref.E3 for detail). 
Table | shows the quasiparticle residues at three different k points and the energy shift due 
to both the hole-spin wave and hole-phonon interactions. a™y c 2 K , 2 \ in the fifth column is 
the quasiparticle residue at the point (7r/2,7r/2) calculated by perturbation theory. The 
numerically computed mass renormalization factor due to the hole-phonon interaction is 
Anum- As expected, the quasiparticle residue decreases as the coupling constant increases. 



At the (0, 0) point the hole loses its quasiparticle property more rapidly, as the vanishingly 
small spectral weight indicates. This is because the state at the (0, 0) point is located at 
a high hole energy region and hence it is more vulnerable to additional decay from the 
hole-phonon interaction. The energy shift due to the hole-phonon interaction is generally 
small (—0.975 for g = 2.0 J), compared with that from the pure hole-spin wave interaction 
(—5.378). This suggests that a perturbative treatment for the hole-phonon interaction is 
possible. According to Eq. \Q the hole self-energy is composed of two terms, one from the 
hole-spin wave interaction and the other from the hole-phonon interaction. Thus the self- 
energy can be written as E(fc,u;) = T^ m (k,u) + S 9 (/c, u), where the former comes from the 
hole-spin wave interaction and the latter from the hole-phonon interaction. If we define X 
and Y as follows 

d 

X = 1- —ReE m (k,uj)\ g=0 , 

d -> -> -> 

Y = --^ Re + £ m (£, u) - £ m (£, u) | g=0 ) , 

the quasiparticle residue can be approximated by 

+ Y X [ 

i, ,y, ,y 



The X, Y are determined by 



x ll-(j) + (j?]- (5) 



A" 



a u (g = 0) 

y = l 1 



H a k(9 = 0) 

Substituting X and Y into Eq. |5| leads to the a("/ c 2>vr /2) i n the fifth column of Table |. For 
coupling constants less than 1.5J, agreement with the numerical results is quite satisfactory. 
However, the perturbative treatment based on the two term expansion breaks down for 
g > 1.5J, since (Y/X) increases significantly for the strong hole-phonon interaction. In fact, 
Y is the approximate hole-phonon mass renormalization constant A evaluated numerically 
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on the basis of the t — J model Hamiltonian. This is listed as \ num in the sixth column. 
The mass renormalization constant for a small coupling constant can also be computed 
from perturbation theory. We used a similar method employed by Ramsak et alR Since 
the hole-phonon coupling constant is small compared with the hole-spin wave interaction 
strength, the effective mass of the hole can be computed using the lowest-order perturbation 
correction to the hole self-energy, 

i r ,9 a k-g9 2 



where uz(g) and are the quasiparticle dispersion functions in the presence of the hole- 
phonon interaction and in its absence, respectively. Approximately, is given by 

(k- (±7r/2,±7r/2)) 2 



2m eff 

where m e ff is the averaged effective mass of the hole near the bottom of the energy dispersion 
function, which is taken as 3.36/t according to Martinez et a/.i. From the second derivative 
around the points k = (±7r/2, ±7r/2), we arrive at 

1 1 ^ lQag 2 m e ff r n ^ q 3 

m ef f(g) m ef f n Jo (q 2 + 2m e// fi) 3 

2™ „ . rv m „,3 



1 2ag z m e ffr. ptm y 



rUm y° 

i 4 i "yj^rw^ < 7 ' 



m ef f vrfi Jo (y 2 + l) 3 
where y m is given by ir / ^2m e //0 and a is the averaged quasiparticle residue near the point 
(7r/2,7r/2). Hence, 

Ae = m eff(g) _ l 1 l 

e m eff 1 — 1.5ag 2 m e ff /(tt£1) 

w 1.5ag 2 m ef f/(iTQ) , (8) 

where a small A e // is assumed. For g/J = 0.5 and 1.0, the substitution of a ~ 0.34 gives 
0.127 and 0.822, respectively. The former value is favorably compared with 0.102, the 
numerically calculated mass renormalization constant for g/J = 0.5. Clearly g/J = 1.0 is 
too strong to apply perturbation theory, because of the too large effective mass, m e ff = 



3.36/t. Besides which, the large value of J2m e ftQ, namely, 1.16 makes the calculation very 



sensitive to the upper bound, y m , for the integration, leading to an additional difficulty. 
A straightforward calculation of X e ff for noninteracting electrons yields 0.016, 0.068, 0.167 
and 0.342 for g/J = 0.5, 1.0, 1.5 and 2.0, respectively. In this case, the above mentioned 
difficulties do not occur because of the small effective mass, wi e // = 1/21 The calculation 
shows that the mass renormalization factor is more enhanced for the strongly correlated 
electrons than the factor for the noninteracting electron system. This enhancement factor 
which is 6.4 for g/J = 0.5 is substantially large, compared with 3.5 which Ramsak et a/.0 
reported. The discrepancy between these two values originates from the use of different 
definitions for the mass enhancement parameter. The authors of this paper obtained the 
enhancement factor by explicitly computing the change in the curvature of the quasiparticle 
dispersion along the line (tt/2,tt/2) — (0,0), namely, An = m\\(g)/m\\ — 1, while the effective 



mass in the present calculation, A e // = y rri\\(g)m^(g) /m\\m± — 1, is averaged out around the 

point (vr/2, vr/2). This indicates even larger mass enhancement along the line (0, tt) — (0, 0), 
which is consistent with the observation that the hole moves slowest along this direction. 

III. INCOHERENT SPECTRUM AND LUTTINGER'S THEOREM 

Because of the additional scattering channel, we expect that there are more spectra as- 
sociated with the optical phonon excitations above the quasiparticle pole in the spectral 
density function. Fig. |2| shows the spectral density function for four different hole-phonon 
coupling constants g. As g increases, the spin wave peaks become more suppressed, whereas 
the phonon peaks become stronger. Since the coupling strength for the hole-spin wave inter- 
action is much larger than that for the hole-phonon interaction, the weak phonon features 
appear on top of the spin wave peaks as a satellite structure. For g > 1.5 J, the phonon 
induced peaks are even sharper and larger in height than the spin wave peaks. Especially 
the appearance of the multiple phonon peaks below the chemical potential for g = 2.0J, 
compared with the spectral density function for g — 0, is clearly noticeable. The spikes due 
to finite size effects are gradually reduced, since those artificial peaks get smeared out due 
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to the additional decay induced by the hole-phonon interaction. Fig. |3| presents the hole 
density of states for various coupling constants. The general behavior for the hole density 
of states is similar to that for the spectral density function, since in most of the Brillouin 
zone the spectral density function is very similar to that at the point (7r/2,7r/2). As the 
hole-phonon coupling constant increases, the position of the strongest phonon peak right 
above the Fermi energy (uo = 0) is shifted upward, although this is not well pronounced due 
to finite size effects. This is understood based on a (lattice) polaronic formation. For the 
strong hole-phonon coupling constant, the hole is surrounded by an increasing number of 
phonons. Hence, approximately 2 phonons are involved in the formation of the polaron for 
g = 2.0 J, since the position of the peak is close to 1.0 J and the optical phonon frequency is 
0.5J. The figure for g = shows, however, 2.5 spin waves with energy 2J (i.e. from the top 
of the spin wave band where the density of spin wave states is sharply peaked) participate 
in a magnetic polaron. The incoherent spectrum below the Fermi energy in the density of 
states is crucial in satisfying Luttinger's theorem0, as far as the quasiparticle (coherent) 
contribution to the density of the occupied states (at T = 0) is substantially suppressed due 
to the strong hole-spin wave (and phonon) coupling. The momentum distribution function 
is shown in Fig. We chose a much larger cluster (240 x 240) for n(k) calculation, than 
the original cluster (24 x 24) used for the self-consistent calculation of the self-energy in 
the k points along the line (n/2, tt/2) — (0, it) or S — Y. This provides detailed information 
about the behavior of n(k) in the vicinity of the Fermi surface. The four very elongated 
ellipses in the inset denote the Fermi surface. The distribution function shows a sharp drop 
at the same k point for all the coupling constants we have studied, as seen in the figure. 
We also compared the doping concentration x computed from the spectral density function, 
with the ratio of the number of k states inside the Fermi surface to the number in the en- 
tire (antiferromagnetic) Brillouin zone. The former yields x = 0.030, 0.031, 0.032, 0.030, 
0.030 for g = 2.0J, 1.5J, 1.0J, 0.5J, 0, respectively, while the latter shows 0.032. They 
agree with each other within less than 4% on the average. These two features numerically 
verify Luttinger's theorem in the t — J — g model for a small doping concentration. As 
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the coupling constant increases, the distribution function inside the Fermi surface decreases 
gradually from 0.36 to 0.23. This is due to the reduced quasiparticle residue for the strong 
hole-phonon interaction, as seen in Table |. At the same time, however, some density of the 
occupied states also appears outside the Fermi surface. This is because in the t — J — g model 
the spectral density function A(k, uj) at the k points outside the Fermi surface possesses a 
strong incoherent tail below the chemical potential. 



IV. VERTEX CORRECTIONS 

Since quasiparticles move coherently on a reduced energy scale 2 J, the Fermi energy is 
quite small for a small doping case, i.e. almost on the order of the phonon frequency or even 
less, which signals a possible breakdown of the standard strong (phonon) coupling theory. 
In the present section, the Migdal-type vertex corrections are studied in the t — J — g model. 
Below k is defined as (k, ik n ) where k n is a Matsubara frequency. Hence a summation over k 
means the summation over both momenta k and Matsubara frequencies k n . The lowest order 
vertex corrections to the hole-phonon interaction in Fig. |5| can be written as gT(k, k + q), 
where 

T(k, k + q) = -^-J2 G(k')G(k' + q)B(k - k') . (9) 

Bik — k') is the Green's function of the optical phonon. Using the spectral representation 
for the hole Green's function and converting the summation over Matsubara frequencies into 
a contour integration leads to 

1 



T(k, k + q)=g 2 J2 / dudu'A{k', uj)A{k' + q, uj') 



uj — UJ' 



x { ., * , [*V) - N(Q) - 1] - T 7~ 7,[F(uj) + N(Q)] 

—ik n + uj + s I —ik n + uj — i I 

1 -\F(uj')-N(n)-l]+ ,, * \F(uj')+N(m, (10) 



— ik n + uj' + Q —ik n + uj' — Q 

where F(u) and N(uj) are the Fermi and Bose distribution functions respectively, and the 
iq n — ► limit is taken first for numerical simplicity. Since the numerical computation for 
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q 7^ shows a similar result to q = case, we restrict ourselves to the latter case. By taking 
the T = and q — > limits and the analytic continuation zfc n — > ko + i8 as well as by noting 
that Q(x) — 1 — 0(— a;) and that (an integral is in the principal value sense) 

,A(%,<J) 



J UJ — uj' 



we arrive at 



r(Jfe„) = 2g 2 J2 J duA(k,u)ReG(k, 



r ©(w) e(-w) , 

x { — 1 i } 

— ko + uj + — id —k + uj — — i5 



2g 2 J2 J duA(k,u)ReG(k,u) 



1 / X 

x — ; 7T~- — r? • (11) 

— Ko + ^ + \ lsign(uj) — id 



Therefore, the real and imaginary parts of the lowest order vertex correction are found as 

dujA(k,uj)ReG{k,uj) 

-co 



x (12) 

— ko + uj + Vtsign{uo) 



and 



ImT(k Q 



2ng 2 Eg A(k, k + Q)ReG(k, k + Q) if k < -Q 
if -Q < ko < Q 

2ng 2 Eg A(k, k - Q)ReG(k, k - Q) if k > Q . 

Fig. |6] and Fig. [7] show the real and imaginary parts of the lowest order vertex correction 
for several hole-phonon coupling constants, respectively. As the coupling strength increases, 
the real and imaginary parts of the vertex correction grow. In spite of the expectation 
that the Migdal approximation may break down due to the small hole band width 2J 
determined from the t — J model0, the vertex correction is much smaller than unity for 
up to g ~ 2 J. For a noninteracting electron system, first order vertex correction has been 
known to be of the order of uup/Ep^, where u)d is the Debye frequency and Ep is the 
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Fermi energy. But, the adiabatic argument valid for a weakly interacting system breaks 
down for strongly correlated electrons, since the Fermi velocity is comparable to or even less 
than the phonon phase velocity in a considerable part of the Brillouin zone. According to 
Eq. first order vertex correction is roughly proportional to the square of the quasiparticle 
residue. This indicates a possibility that the significant renormalization (0.2 — 0.3) of the 
quasiparticle residue for a strongly correlated electron system makes the vertex correction 
much reduced, thereby accounting for the small vertex correction from the hole-phonon 
interaction in the t — J — g model. Hence, the present calculation numerically corroborates 
using the noncrossing approximation for the self-energy both in the hole-magnon and hole- 
phonon interactions. 

V. OPTICAL CONDUCTIVITY IN THE t-J-g MODEL 

The current operator after Bogoliubov transformations of the spin wave operators 
becomes^ 

2et 

JM = -jf £ hh}[C %>p al p ^+ D^_ n+f+q l , (13) 

k,p 

where the bare current vertices Cg^,, D^p are defined as 

C k,p = U k~p Sin P* + V k-p Sili k * » 

D k,p = v k-p sin Vx + sin k x . (14) 

Due to the special nature of the interaction vertices, it was established that the lowest order 
contribution to the optical conductivity dominates at uj ^ 0c!. Hence, we consider only the 
lowest order diagrams in the present study. The lowest order contribution comes from two 
diagrams in Fig. |8| owing to the structure of the current operator Eq. |T^. In the combined 
limits of zero temperature (T — > 0) and long wavelength electromagnetic radiation (q — > 0), 
the optical conductivity becomes 

71 2t fi° -» 

0i (9o) = — (^rr) 2 / dw y ^A{k 1 u) - g ) 
g A Jo ^ 

13 



x E^b>-^). (15) 

p 

Fig. ^| presents the optical conductivity for the t — J — g model, as the hole-phonon 
coupling constant varies. As expected from the corresponding spectral density function, 
the contribution to the conductivity from the multi-magnon excitations ("string structure") 
decreases, while strong absorption appears right above the 2J peak. This new peak in the 
absorption comes from the hole-phonon interaction, as can be seen in the corresponding 
spectral density function. Generally, the peak at 2J peak and higher energy incoherent 
spin wave peaks are suppressed and broadened in the presence of the strong hole-phonon 
interaction. This may be associated with the growth of featureless spectral weight at the 
mid-infrared region, when the Cu02 plane is doped with charge carriers^. 



VI. CONCLUSION 

The influence of the hole-phonon interaction on various physical quantities was studied 
within the noncrossing approximation for the spin wave and optical phonon interactions 
on the same footing. As the hole-phonon coupling constant g increases, the quasiparticle 
residue is further reduced and spin wave peaks in the spectral density function and optical 
conductivity are more suppressed. Phonon peaks in the spectral density function A(k,u), 
instead, grow more pronounced around the quasiparticle pole at a low energy on the scale of 
Q . A sharp drop in the hole momentum distribution function is found for all the hole-phonon 
coupling constants we have studied. The invariance of the volume enclosed by the Fermi 
surface for all the chosen hole-phonon coupling constants numerically verifies Luttinger's 
theorem for doped holes in the t — J — g model. Our numerical estimate of the lowest order 
vertex corrections to the hole-phonon coupling vertex g due to the hole-phonon interaction, 
gives relatively small values <C 1. This means that electron-phonon vertex corrections are 
not important and Migdal's approximation can be used in calculating the hole self-energy. 
The smallness of the effective hole bandwidth of order 2J, is compensated by suppressed 
quasiparticle residues. Due to the presence of additional phonon induced absorption, 2 J and 
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higher hole-multi-spin wave peaks in optical conductivity are suppressed and broadened. 
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TABLES 



g/J 


°(o,o) 


a (0,7r/2) 


a(n/2,n/2) 


a (*/2,w/2) 


A^/J 


Anum 


0.0 


0.103 


0.408 


0.349 


0.349 


-5.378 


0.000 


0.5 


0.001 


0.376 


0.337 


0.337 


-5.425 


0.102 


1.0 


0.002 


0.322 


0.304 


0.305 


-5.603 


0.424 


1.5 


0.003 


0.258 


0.266 


0.274 


-5.934 


0.894 


2.0 


0.007 


0.231 


0.216 


0.266 


-6.353 


1.765 



TABLE I. Quasiparticle residue, energy shift and mass renormalization factor due to the 
hole-phonon interaction. The first three columns after the column corresponing to the value of 
g/J are quasiparticle residues at three different k points and the fourth column is the energy shift 
due to the hole-magnon and hole-phonon interactions at the (ir/2,ir/2) point. a("/ c 27r / 2 ) is the 
quasiparticle pole strength at the (7r/2, it/2) point obtained by perturbation theory. \ num is the 
numerically calculated renormalization factor due to the hole-phonon interaction. 
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FIGURES 

FIG. 1. Four different diagrams contributing to the self-energy, magnon emitting and absorbing 
processes, and phonon emitting and absorbing ones. 

FIG. 2. Spectral density function for four different coupling constants g at k = (7r/2,7r/2) for 
a 24 x 24 momentum cluster. The figures for g ^ are shifted upward by 0.2. The doping 
concentration is 3.2%. The phonon frequency Q is equal to 0.5 J. 

FIG. 3. Density of states for the dressed holes for four different coupling constants g for a 24 x 
24 momentum cluster. The figures for g ^ are shifted upward by 0.2. The doping concentration 
is 3.2%. 

FIG. 4. Momentum distribution function for four different coupling constants g for a 24 x 24 
momentum cluster. It is scanned from S= (ir/2,ir/2) to M= (vr,0). The ellipse-like Fermi surface 
pockets are calculated numerically for the doping concentration 3.2%. 

FIG. 5. Lowest order vertex correction due to hole-phonon interaction. 

FIG. 6. Real part of the lowest order vertex correction due to hole-phonon interaction for a 
24 x 24 momentum cluster. The figures for g/J > 1.0 are shifted upward by 0.1. The doping 
concentration is 3.2%. 

FIG. 7. Imaginary part of the lowest order vertex correction due to hole-phonon interaction for 
a 24 x 24 momentum cluster. The figures for g/J> 1.0 are shifted upward by 0.1. The doping 
concentration is 3.2%. 

FIG. 8. (a) Lowest order diagram for the optical conductivity, corresponding to an intermediate 
magnon emitting process, (b) Other lowest order diagram for the optical conductivity. At zero 
temperature, this contribution vanishes, since it involves a magnon absorption. 
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FIG. 9. Optical conductivity calculated by the analytic expressions given in the text for four dif- 
ferent hole-phonon coupling constants g for a 24 x 24 momentum cluster. The units are 27re 2 x 10 2 
choosen for convinience (compare Ref. 9). The figures for g ^ are shifted upward by 20. The 
doping concentration is 3.2%. 
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